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1. INTRODUCTION 

c 

The Burst and Transient Source Experiment (BATSE) is one of 
four instruments on the Compton Observatory which was launched by 
the space shuttle Atlantis on April 5, 1991. As of mid-March, 
1994, BATSE detected more than 925 cosmic gamma-ray bursts and 
more than 725 solar flares. Pulsed gamma-rays have been detected 
from at least 16 sources and emission from at least 28 sources 
(including most of the pulsed sources) has been detected by the 
Earth occultation technique. UAH participation in BATSE is 
extensive but can be divided into two main areas, operations and 
data analysis. The daily BATSE operations tasks represent a 
substantial level of effort and involve a large team composed of 
MSFC personnel as well as contractors such as UAH. The 
scientific data reduction and analysis of BATSE data is also a 
substantial level of effort in which UAH personnel have made 
significant contributions. 

2. MISSION OPERATIONS 

W. Paciesas served as BATSE Mission Operations Software 
(MOPS) Development Manager, chairing the Level V Configuration 
Control Board for the MOPS software. Paciesas also supervised 
distribution of the Individual Burst Data Base (IBDB) files to 
co-investigators at Goddard Space Flight Center (GSFC) and the 
University of California, San Diego (UCSD) , and to the CGRO 
Science Support Center. K. Squier provided maintenance of 
software to produce summary reports of Individual Burst Data Base 
(IBDB) contents for selected triggers. Squier also continued 
development of software to update certain IBDB contents and/or 
produce additional IBDB files without regenerating the entire 
IBDB for a particular trigger. Z. Shariff and G. Richardson 
assisted in BATSE data operations on a regular basis. Shariff, 
Richardson, T. KosfVut, R. Mallozzi, and J. Brainerd performed 
burst operations regularly. 

Paciesas led an effort to plan replacement of the MOPS burst 
catalog with an improved trigger catalog. A new structure and 
organization for the catalog was developed in consultation with 
various BATSE science team members. Paciesas wrote a Software 
Change Plan to implement the catalog with a minimum of impact on 
normal operations. Squier devloped code for initial conversion 
of the old MOPS catalog into the .new format. 
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V. Chaganti assumed responsibility for maintenance of the 
HER/SHER correction software. His first task involved correction 
of the first HERB/SHERB readout records and merging of these into 
the corrected HER/SHER files. 

4. BURST DATA ANALYSIS 

Paciesas, Pendleton, Mallozzi and Koshut were part of the 
collaborative effort to produce the first BATSE catalog of gamma- 
ray bursts. 1 Pendleton was primarily responsible for calculation 
of burst fluxes and fluences, while Koshut had primary 
responsibility for calculation of burst durations. The durations 
were shown to have a bimodal distribution 2 with a minimum around 
1—2 s. To investigate this further, Koshut led a preliminary 
study contrasting the properties of short events with short 
spikes in long events. 3 

Pendleton, Koshut and Chaganti developed software to be used 
for production of the second burst catalog. Improvements to the 
existing algorithms for calculation of burst fluxes, fluences, 
and durations were implemented and a user-friendly interface 
developed. The ultimate goal is produce these numbers routinely 
as part of daily burst operations rather than as a later 
production data analysis task. Chaganti and Shariff completed 
production of fluxes and fluences for the second catalog using 
this software. Koshut and Mallozzi completed production of most 
time duration measurments using this software. In a related 
effort for the second catalog, Koshut led a study of systematic 
errors in the location of bursts using the MAXBC datatype. 4 

Paciesas coordinated BATSE spectral analysis efforts among 
UAH, MSFC , GSFC, and UCSD and supervised the time-resolved burst 
spectroscopy key project. Visual searches for spectral features 
in BATSE bursts were continued as part of this project. No 
convincing evidence has yet been found for such features; 5 a 
candidate line found in one burst is probably spurious. 6 ' 7 The 
lack of line features cannot yet be considered inconsistent with 
other observations. 8 ' 9 An automated procedure for exhaustive 
searches for line features is being implemented. 10 Studies of 
burst continuum spectra showed differences among bursts as to 
their temporal evolution. 11-14 The spectroscopy collaboration 
also produced a catalog of burst spectra. 15 
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Pendleton led studies of continuum spectral characteristics 
of gamma-ray bursts using low-energy-resolution data (4 
channels) 16 and medium-energy-resolution (16 channels). 17 
Pendleton collaborated with GRO Guest Investigator (GI) E. 
Fenimore (Los Alamos Nat. Lab.) and others on deriving and 
interpreting the number/ intensity distribution of gamma-ray 
bursts obtained by combining BATSE and PVO data. 18 

UAH personnel were involved in various studies of burst sky 
distributions. Brainerd led a collaboration which found no 
significant clustering of bursts using a nearest-neighbor 
analysis. 19 Koshut and Mallozzi worked with other BATSE team 
members to put limits on heliocentric models for bursts. 20 
Pendleton and Paciesas collaborated with J. Hakkila (Mankato St. 
U. ) and others on constraining galactic burst models from BATSE 
data. 21 

Pendleton, Paciesas and Koshut collaborated with C. 
Kouveliotou (MSFC/USRA) and others in analysis of recurrent 
events detected by BATSE from the soft gamma repeater 
SGR 1900+14. 22 Pendleton and Paciesas also collaborated in 
studies led by Kouveliotou of events from SGR 1806-20 23 and the 
unusually intense event GRB930131. 24 

Paciesas worked with I. Mitrofanov (Space Research Inst., 
Moscow) , J. Norris (NASA/GSFC) and others on searches for time 
dilation effects in bursts which would be expected if they are 
cosmological in origin. The results so far are suggestive but 
not conclusive. 25,26 

Brainerd developed a model for burst spectra from sources at 
comological distances based on Compton scattering in a relatively 
thick obscuring medium 27 ' 28 and studied effects which might mimic 

• • • • , . ■ OQ 

cosmological time dilation of bursts in Euclidian space. * 

5. NON-BURST DATA ANALYSIS 

Paciesas and Mallozzi worked on occultation analysis of 
selected sources. Mallozzi searched for low-energy gamma-ray 
emission from a set of active galactic nuclei (AGNs) detected by 
EGRET. 30 Paciesas led analysis of temporal and spectral 
variability of the quasar 3C 2 73 31 and the transient source Nova 
Muscae. 32 Paciesas analyzed the intensity dependence of Cygnus 
X-l spectral states around a period of unusually low intensity in 
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January-February 1994. 33 Paciesas and Pendleton collaborated 
with various other BATSE team members on studies of the galactic 
black-hole candidate GX 339-4, the transient sources GRO J1719-24 
and GRS 1009-45, and sources in the galactic center region. 34-36 
Paciesas collaborated with the effort led by S. N. Zhang (USRA) 
to develop a method for producing images from occultation 
data. 37-39 Paciesas collaborated with J. Greiner (MPE Garching) 
on analysis of the transient GRS 1915+105 40 and with J. Ling 
( JPL) on analysis of Cygnus X-l. 41 

M. Stollberg continued analysis of data from the x-ray 
binary EXO 2030+375 as part of his Ph. D. dissertation. 42 
Stollberg also assisted the BATSE team in production pulsar 
analysis and worked on documentation of pulsar timing analysis 
algorithms. Stollberg and Pendleton collaborated with R. Wilson 
(MSFC) and others in analysis of the transient pulsar GRO J1008- 
57. 43/44 Pendleton collaborated with Wilson and others on a 
study of intensity/torque correlations in Her X-l. 45,46 Paciesas 
collaborated with B. Rubin (USRA) and others on an investigation 
of the long-period pulsar 4U 1538-52. 47 

Paciesas and Pendleton continued to work with various other 
CGRO Guest Investigators in ongoing studies. These included A. 
Parsons (NASA/GSFC) , R. Edelson (Univ. of Iowa) and M. Machado 
(UAH) . 

5 . OTHER ACTIVITIES 

Brainerd and Pendleton served on the local organizing 
committee for the Gamma-Ray Burst Workshop held in Huntsville on 
20—22 October 1993. Brainerd co-edited the workshop proceedings, 
along with G. Fishman (MSFC) and K. Hurley (UC Berkeley) . 
Paciesas, Pendleton and Mallozzi presented contributed papers at 
the Second Compton Symposium in College Park, MD, during 20—22 
September 1993. Paciesas presented an invited talk on 
"CGRO/ BATSE occultation studies of galactic black hole candidates 
and AGNs" at the workshop on Pairs , Gamma-Rays and Black Holes in 
Koninki, Poland, during 5-8 October 1993. Stollberg and Paciesas 
presented contributed papers at the conference on The Evolution 
of X-Ray Binaries in College Park, MD, during 11-13 October 1993. 
Paciesas presented invited lectures on "BATSE Observations of 
Extragalactic Objects" and "BATSE Observations of Galactic 
Compact Objects" at the NATO Advanced Study Institute on The 
Gamma-Ray Sky from Compton /GRO and Sigma in Les Houches, France, 
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during 25 January-4 February 1994. 


Paciesas continued to serve as BATSE representative to the 
CGRO User's Committee, attending the meeting on 23-24 September 
1993. Paciesas reviewed and provided updated BATSE inputs to the 
CGRO Project Data Management Plan and to the NASA Research 
Announcement for phase 3 of the GI program. 

Meetings of the Burst Spectroscopy Team were held on 9—10 
September 1993 at GSFC and on 17-18 February 1994 at UCSD (both 
were attended by Paciesas and Pendleton) . 

Copies of selected publications involving UAH personnel as 
principal author are attached. 
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THE NEAREST NEIGHBORS OF BATSE GAMMA-RAY BURSTS: 
NARROWING THE POSSIBILITIES 
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ABSTRACT 

The large errors in the locations of gamma-ray bursts observed by the 
Bursts and Transient Source Experiment (BATSE) forces one to adopt statis- 
tical comparisons of the burst catalog to models of burst isotropy to determine 
if bursts are repeating or clustered. In the first BATSE catalog 1 , a nearest 
neighbor analysis finds a deviation from isotropy. 2 In a recent article 3 It was 
shown that known instrumental effects cannot produce significant small scale 
anisotropies. It was also shown that burst repeater models can produce the 
observe anisotropies. In this paper we examine in more detail repeater burst 
models. We also show that nearest and farthest neighbor analyses of more re- 
cent bursts fail to find significant small scale clustering. 

INTRODUCTION 

A nearest neighbor analysis 2 of the first BATSE Gamma-Ray Burst Cata- 
log 1 finds a maximum deviation from the isotropic nearest neighbor cumulative 
distribution of 2% statistical significance. 3 This contradicts an earlier analysis 
that examines the average value of the nearest neighbor separation and finds 
no significant deviation from isotropy. 4 The question arises whether the result 
of the first named analysis is a consequence of statistical fluctuations or of ei- 
ther instrumental or physical processes. The authors of one article 5 assert that 
the deviation must be of instrumental origin, because a similar farthest neigh- 
bor analysis of the first Catalog finds a maximum deviation from the isotropic 
cumulative distribution of 11% significance, which they regard as statistically 
significant and unphysical. 

In a recent article 3 it was shown that there are no systematic effects that 
can produce a deviation of the magnitude seen in the data. The effect of the 
blockage of the sky by the earth produces a maximum deviation of D = 0.006 
from the isotropic cumulative distribution, which is much smaller than the 
D = 0.12 deviation of the First Catalog. While one might expect a small 
angle anisotropy from the flux dependence of the sky map to introduce a devi- 
ation, we find from Monte Carlo simulation that the maximum deviation from 
the isotropic cumulative distribution is D = 0.01 without any spacecraft re- 
orientations and < 10 -3 with the reorientations experienced by the Compton 
Gamma-Ray Observatory (CGRO). Moreover, performing the nearest and far- 
thest neighbor analyses in the spacecraft coordinate frame produces maximum 
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deviations of « 50% significance. If the deviation in the First Catalog is non- 
statistical, it must have a physical origin. 

The suggestion of Lamb and Quashnock 2 that repeating burst sources are 
responsible for the nearest neighbor results is not invalidated by the farthest 
neighbor analysis. 3 In fact, if repeating burst sources are confined to the galac- 
tic plane, the resulting nearest and farthest neighbor distributions agree with 
the First Catalog distributions without violating the limits on the dipole and 
quadrupole moments. 

In this paper we discuss disk confined repeater models in more detail, 
presenting their nearest and farthest neighbor cumulative distributions. We also 
present nearest and farthest neighbor analyses of a recent set of 260 gamma-ray 
bursts; the results are consistent with an isotropic distribution of bursts. 


BURST MODELS 


The nearest neighbor of a given burst is primarily determined by the local 
density of bursts on the sky. The magnitude of the density fluctuations therefore 
determines the shape of the cumulative distribution. On the other hand, both 
the magnitude and topology of density fluctuations strongly affect the farthest 
neighbor cumulative distribution. This implies that the farthest neighbor distri- 
bution is much more model dependent than the nearest neighbor distribution. 

Two simple models illustrates this point. First, imagine that aN bursts 
are distributed on the upper celestial hemisphere and (1 — a) N bursts are dis- 
tributed in the lower celestial hemisphere, where N, the total number of bursts, 
is very large. The nearest neighbor distribution is then 


( 1 — cos#\ . ( 1 — cosfl \ 

f " =aeXP (-^^J +(1 ~ a)eXP U 4iY(l- 0 ) J 


and the farthest neighbor distribution is 

_ /, \ / 1 — cos 9 

F f = (l-a)exp I ) + a exp 


1 — cos 9 
4JV (1 — a) 


( 1 ) 

( 2 ) 


The isotropic nearest and farthest cumulative distributions (a = 1/2) equal 
1/2 when 1 — cos 9 = 2iVln2. For this value, F n > 1/2 and Ff < 1/2 for all 
a ^ 1/2. Now consider a second model in which aN / 2 bursts axe distributed on 
one half of the upper hemisphere and the same number of bursts are distributed 
on the geometrically opposed half of the lower hemisphere. The remaining two 
half-hemispheres have ( 1 — a) N/2 bursts apiece. In this model the nearest and 
farthest cumulative distributions are both given by equation (1), proving the 
importance of the topology of density variations. 

These simple models illustrate a second point: the nearest neighbor cumu- 
lative distribution lies above the isotropic distribution for small vales of 1 — cos 9 , 
but the farthest neighbor distribution can lie on either side of the isotropic dis- 
tribution. However, the magnitude of the deviation from isotropy of the nearest 
neighbor distribution also depends on the density gradient, for a must be less 
than « 0.15 to produce a deviation as large as in the left hand plot of Figure 1. 
This effect disappears as N increases. 

The source of the anisotropy also determines whether the deviation of the 
cumulative distribution from the isotropy distribution changes as the number 




Figure 1. The nearest and farthest neighbor cumulative distribu- 
tions: effect of sample size. Theoretical curves for a model in which 
burst sources are isotropically distributed and a fraction of the burst 
sources producing 4 bursts apiece. The isotropic cumulative distri- 
bution is given by the dotted line, the model nearest neighbor distri- 
bution by the solid line, and the model farthest neighbor distribution 
by the dashed line. Left: 200 single burst sources and 15 multiple 
burst sources. Right: 400 single burst sources and 30 multiple burst 
sources. 


of bursts increases. If the anisotropy arises from the presence of a disk compo- 
nent, then the anisotropy persists. If repeating sources produce the anisotropy 
through local density enhancements over a location error angular scale, then the 
cumulative distributions should go to the isotropic distribution as the number of 
bursts increases beyond the point where the location errors about each repeating 
source overlaps with adjacent repeating sources. Figure 1 demonstrates the de- 
pendence of the cumulative distribution on the sample size. As a consequence, a 
statistical search for repeaters should use consecutive sets of bursts with a fixed 
number of bursts in each set. 

Figure 1 shows that isotropically distributed repeating sources produce a 
farthest neighbor cumulative distribution that falls below the isotropic distri- 
bution. If the repeaters are confined to the galactic plane, as in Figure 2, the 
farthest neighbor cumulative distribution can fall above the isotropic distribu- 
tion. 

The model used in Figure 2 has two components, with one component an 
isotropically distributed set of single burst sources and the second component 
a set of multiply bursting sources confined to the galactic plane. Such models 
can reproduce the nearest neighbor distribution of the First BATSE Catalog 





Figure 2. The nearest and farthest neighbor cumulative distribu- 
tions: quadruple repeaters versus double repeaters. Two hundred 
single burst sources are isotropically distributed. Additional repeat- 
ing sources are uniformly distributed between ±10° galactic latitude. 
The position of each outburst relative to the source is given by a 
gaussian distribution function with average offset of 5°. The lines 
are defined as in Fig. 1. Left: 15 repeating sources burst 4 times 
apiece. Right: 30 repeating sources burst 2 times apiece. 


without violating constraints on the farthest neighbor distribution or the value 
of (sin 2 6). The right hand plot of Figure 2 deviates from the the observed 
nearest and farthest neighbor distributions with significances of 68% and 31% 
respectively. In this model one finds sin 2 b — 0.26, which is less than 3 a from 
the value in the First Catalog. 

RECENT RESULTS AND CONCLUSIONS 

The Kolmogorov- Smirnov statistic for the nearest neighbor distribution 
from the First BATSE Catalog 3 of 260 bursts is N X / 2 D = 1.86, which has a 
2% significance. A farthest neighbor analysis gives N l ! 2 D = 1.56, which has an 
11% significance. 5 These values are consistent with both a statistical origin and 
a physical origin. 

A more recent ensemble of 260 gamma-ray bursts disjoint from the first 
catalog shows no evidence for anisotropy. Their nearest and farthest neighbor 
cumulative distributions are given in figure 3. The maximum deviation of the 
nearest neighbor cumulative distribution is below the isotropic cumulative dis- 




Figure 3. The nearest and farthest neighbor cumulative distribu- 
tions: recent gamma-ray burst sample. A set of 260 gamma-ray 
bursts is selective by dropping all gamma-ray bursts with locations 
by the BATSE data type MAX-BC from the set of triggers between 
1819 and 2494 inclusive (GRB920814b to GRB930817). These bursts 
will be part of the 2B and 3B gamma-ray burst catalogs. The lines 
are defined as in Fig. 1. 


tribution, and gives N l / 2 D = 0.77, with significance 83%. The farthest neighbor 
distribution has N l / 2 D = 0.85, with significance 75%. Again the maximum de- 
viation is below the isotropic cumulative distribution. This provides no support 
for the belief that the deviation found in the first catalog has a physical origin. 
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THE COMPTON ATTENUATION MODEL 
OF COSMOLOGICAL GAMMA-RAY BURSTS 
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ABSTRACT 

A gamma-ray burst spectrum can be modeled by passing a power law 
spectrum through a medium that is optically thick to Compton scattering. The 
roll-over of the Klein-Nishina cross section at high energies produces a break 
in the spectrum at several hundred keV if the gamma-ray source is at z ss 1. 
Photon-photon pair production limits the scattering region’s size to greater than 
« 0.1 parsec. Because the optical depth is greater than unity for a length scale 
of several parsecs, the density must be of order 10 5 cm -3 , which occurs only 
in molecular clouds at the centers of galaxies. This model therefore precludes 
source models employing objects that are common to the galactic plane — for 
example, merging neutron stars. The model spectrum is described by five free 
parameters: a cosmological red shift, a Thomson optical depth, a power law 
index, a metallicity, and an amplitude. The attenuation model can be tested by 
comparing the consistency of model parameters derived at different times over 
the duration of a burst. Small angle scattering of x-rays by dust within the 
molecular cloud produces an x-ray afterglow. A consequence of this model is 
that any optical or ultraviolet radiation is heavily absorbed, making detection 
of gamma-ray bursts at these energies unlikely. 

INTRODUCTION 

A long standing problem in constructing models of gamma-ray bursts is 
modeling the non-thermal burst spectrum. These spectra are generally a power 
law above IMeV, and emit their peak power at several hundred keV. 1,2 Be- 
fore BATSE, when the predominate models were emission from neutron stars, 
the small fraction of energy in the x-ray band caused considerable theoretical 
difficulty, and became known as the x-ray paucity problem. Now we have cos- 
mological models which have no predefined characteristic energy, implying that 
the characteristic spectral energy peak is still a problem. 

The root of the difficulty with cosmological models is the large Lorentz 
factors of 7 > 100 required to overcome photon-photon pair creation. This makes 
the characteristic energy in the rest frame of the emitter ill-defined, because the 
value of the Lorentz factor is not physically constrained. The likely radiative 
processes — photoelectric absorption, cyclotron emission, Compton scattering — 
are further unconstrained in optical depth, electron distribution, magnetic field 
strength, and seed photon density and energy. 

These considerations suggest that the characteristic energy of gamma-ray 
bursts is set by processes in the rest frame of the host galaxy and is associated 
with the electron rest mass energy. This article examines Compton attenuation 
as the mechanism responsible for producing the observed spectra. More detailed 
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analyses of this model axe presented elsewhere. 3,4 

THE COMPTON ATTENUATION MODEL 

Because the purpose of this model is to explain the characteristic peak 
values of the vF v curve independent of the Lorentz factor of the source, I must 
assume that the source spectrum is a power law. In the model the source is 
embedded in a cool, static medium with Thomson optical depth much greater 
than unity. Compton scattered radiation can annihilate with other burst pho- 
tons and produce electron-positron pairs. If significant scattering occurs closer 
than ~ d a = 0.04 (r/10) exp (r/2 — 5) pc to the source, where r is the Thomson 
optical depth of the scattering region, the number of pairs created raises the 
optical depth by at least 1, which produces a rapid and permanent extinction 
of the gamma-ray burst. As this contradicts the observations, the scattering 
region must be larger than d s . The travel time between scatterings is therefore 
> 10 6 s, so the scattered radiation does not contribute to the burst. 

A second process affects the x-ray spectrum: photoelectric absorption. 
Because the scattering region is primarily composed of hydrogen, this process 
affects the spectrum at the lower energy limit of the BATSE instrument (ap- 
proximately 20keV). The energy of unity optical depth is pushed far below the 
instrument threshold for a cosmological red shift of z « 1. As a consequence, 
photoelectric absorption is generally unimportant for fitting the spectra of bursts 
observed by BATSE. 

Figure 1 shows the Compton attenuated spectrum for several values of 
the Thomson optical depth. At high energies the scattering optical depth falls 
to zero because of the v~ x dependence of the Klein-Nishina cross section above 
IMeV. The transition from the Thomson cross section to the Klein-Nishina cross 
section produces a “broken” spectrum below IMeV. Below lOkeV the effects of 
photoelectric absorption appear. In the absence of photoelectric absorption the 
spectrum follows a power law with the index of the source power law spectrum. 

The high Thomson optical depth in this model limits it to sources at the 
centers of galaxies. Because the scattering region is large, it must be associated 
with the interstellar medium. But the optical depth through the plane of nearby 
spiral galaxies is much less than unity. This leaves molecular clouds with den- 
sities of order 10 5 cm -3 , which exist only at the centers of galaxies. This model 
is therefore inconsistent with sources common to the planes of galaxies, such as 
merging binary stars. The source must in fact be unique to the cores of galaxies, 
which suggests that the energy source in this model is a massive black hole. 

X-RAY AFTERGLOW 

A small number of bursts have x-ray emission lasting up to several hun- 
dred seconds after the gamma-ray emission falls below background. 5,6,7 The 
afterglow spectrum falls rapidly and can be fit by a black body continuum with 
a temperature of approximately IkeV to 2keV. To this point the only interpre- 
tation offered for this radiation is thermal emission from the burst source. For 
extragalactic sources this implies a source size of « 10 10 cm for 7 = 1, which is 
larger than the size derived from the rise time. 

Because the Compton attenuation model requires a dense molecular cloud 
for the scattering region, it provides a natural explanation for the x-ray after- 
glow — coherent scattering with dust. 4,8,9,10 X-rays produced in the gamma-ray 
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Figure 1. Spectra produced by Compton attenuation of a power law 
with index 8 = —3. The scattering region has 0.5 of the solar abun- 
dance of metals. The Thomson optical depths from top to bottom 
are r = 0, 1, 5, 10, 20, and 30. No cosmological red shift is present. 


burst axe scattered by a small angle into the observer’s line of sight. The longer 
propagation path introduces a delay in the x-ray emission relative to the gamma- 
ray emission. At several keV the scattering angle is of order 10 -4 radian for a 
dust particle diameter of 0.01//m, which produces a time delay of 100 seconds if 
the scattering occurs lOpc from the x-ray source. 

I assume the distribution function for dust grain size is 


/( a ) 


3(4 - m) p m /ao\ m 
47 t a^pg V a / 
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where a < ao is the grain radius, p g is the mass density within each grain, 
and p m is the mass density of metals in the molecular cloud. The distribution 
function is zero when a > Go- The total cross section for coherent scattering of 
x-rays by dust is then 4 
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In this equation the function in brackets is of order unity. For the molecular 
clouds under consideration in this model, one expects the energy of unity optical 
depth to be near lOkeV. For these equations and the differential scattering cross 
section, I ran a Monte Carlo simulation of 10 5 photons to produce Figure 2. 
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Figure 2. X-Ray afterglow. The dotted curve gives the Compton 
attenuated burst spectrum in the absence of dust for r = 10. The 
remaining parameters are those for Fig. 1. The uppermost solid 
curve gives the direct spectrum in the presence of dust, where all 
metals are assumed to be in dust. The dust distribution is given 
in the text. The remaining two solid curves give the delayed x-ray 
spectrum. The upper curve is the delayed emission integrated over 
time 100 s after the burst. The lower curve is the emission integrated 
between 100 s and 10 3 s. 


The dust decreases the direct x-ray flux by a significant factor below 30keV. 
The afterglow spectrum is quite steep above unity optical depth, having a spec- 
trum given by the burst spectrum times E 7 k ~™. At low energies the optical depth 
to scattering exceeds unity, causing the delayed spectrum to be flat. While the 
spectra in Figure 2 resemble thermal spectra with temperatures of « 5keV, 
which is higher than observed, dust scattered spectra at a cosmological red shift 
of « 2 will have the correct shape to fit the observations. One characteristic of 
this model is that the peak of the x-ray afterglow should shift down in energy 
with time. 


OBSERVATIONAL TESTS 

The Compton attenuation theory must satisfy a stringent test: model fits 
to a burst’s spectrum at different times over the duration of the burst must yield 
the same cosmological red shift z. One does not generally expect other models 
to produce a constant value of z. For instance, if the energy of the spectral 
break is dependent on the Lorentz factor 7 of the emitter, and if 7 decreases 




with time, the spectral break energy should fall with time, leading to a z that 
rises with time. The constancy of z with time can also be tested by fitting the 
color-color diagrams of individual bursts 11 with model color-color curves found 
by allowing the power law index to vary. 

A second test of this model is the comparison of the values of z derived 
from spectral fits with the allowed values of z implied by each burst’s position 
on the log iV-log P max curve. High values of z must correspond to low values of 

P max- 

The effects of dust offer additional opportunities for testing this theory. 
The spectrum and time profile of the x-ray afterglow is dependent to some 
extent on the spectrum and time profile of the direct gamma-ray burst. One 
should be able to model this correlation and test for its presence. A second 
test is the presence of an optical flash at the time of the gamma-ray burst. 
Optical emission during the gamma-ray burst will be heavily absorbed by the 
dust, making the observation of gamma-ray bursts at these wavelengths unlikely. 
The observation of an optical flash would pose great physical difficulties for the 
Compton attenuation model. 
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A STUDY OF CONTINUUM SPECTRA OF SHORT-DURATION 
GAMMA-RAY BURSTS OBSERVED BY BATSE 


T. M. Koshut, G. N. Pendleton, R. S. Mallozzi, W. S. Paciesas, M. S. Briggs 
Dept, of Physics, University of Alabama in Huntsville, Huntsville, AL 35899 

ABSTRACT 

We use 4 and 16 channel data from the BATSE Large Area Detectors 
to compare the continuum spectral characteristics of the short- duration bursts 
with those of short spikes within long-duration bursts. We compare frequency 
distributions of the resulting best-fit model parameters for the short bursts and 
for the short spikes within the long bursts. These distributions are different for 
the two types of events. The best-fit model parameters are used to calculate the 
peak energy E ma . x of the v F„ spectrum for each type of event. The resulting 
distributions of E ma , x are consistent within a cosmological scenario if the redshift 
of the short spikes within long-duration bursts is larger than the redshift of the 
short bursts. 


INTRODUCTION 

The existence of two subclasses of classical gamma-ray bursts, defined 
by durations <2.5 seconds and durations >2.5 seconds (see reference 1, and 
references therein), has been confirmed. 2 The time histories of the long-duration 
events often contain substructures with pulse widths that are very similar to the 
durations of the short-duration bursts. These time histories can be very diverse, 
ranging from those that are highly variable over the duration of the burst to 
those having an impulsive spike riding on top of gradually-varying emission. 
By treating each such substructure as a separate event, we hope to investigate 
the possibility that each substructure is a manifestation of the same emission 
mechanism producing the short-duration bursts. 

METHODOLOGY 

BATSE consists of eight uncollimated Large Area Detectors (LADs), span- 
ning an energy range of ~ 20-2000 keV. Details of BATSE instrumentation can 
be found elsewhere. 3 Five BATSE data-types are used in this study: Medium 
Energy Resolution (MERl, Discriminator Science (DISCSC), Preburst (PREB), 
Continuous (CONT), ana LAD Discriminator (DISCLA). MER data consist of 
count rates in 16 energy channels, with a time resolution of 16 ms for the first ~ 
32 seconds after the burst trigger time t 0 , changing to a time resolution of 64 ms 
for the next ~ 128 seconds. MER data begins to accumulate ~ 30 ms after t Q , 
usually making it unsuitable to study the short-duration bursts (which axe often 
over by t 0 ). The PREB and DISCSC data consist of count rates in 4 energy 
channels, with a time resolution of 64 ms, starting ~ 2 seconds prior to t 0 and 
lasting until ~ 240 seconds after t 0 . CONT data consist of count rates in 16 en- 
ergy channels, with a time resolution of 2.048 seconds. DISCLA data consist of 
count rates in 4 energy channels with a time resolution of 1.024 seconds. CONT 
and DISCLA data are only used for the purpose of background subtraction in 
this study. Only the 16-channel data were corrected for dead-time effects; the 



4-channel data would not have suffered any significant dead-time during the 
bursts used in this study. 

A total of 29 short-duration bursts and 13 long-duration bursts were ex- 
amined in this study. A count spectrum was integrated over the duration of each 
short-duration burst. Additionally, within each long-duration burst, a spectrum 
was integrated over short spikes with a pulse width less than 2.5 seconds. This 
criterion resulted in 55 short spikes among the long- duration events. 

The count spectra were deconvolved to produce photon spectra using a 
forward-folding model-fitting technique. A detector response matrix was used 
for each burst, correcting for the detector’s energy response function, angular re- 
sponse function, and for flux scattered by the Earth’s atmosphere and the CGRO 
spacecraft. 4 The photon spectrum model employed in the fit was a Gaussian as 
a function of log 10 (E), given by 


dN A 1 

dE ~ exp ( 2 


( 2 ) 

Three parameters were allowed to vary in this model: A, E ce nt, and cr g . A 
represents the amplitude of the Gaussian, E C ent is the center of the Gaussian, and 
<7 g measures its full-width at half-maximum. Though this model is unphysical, 
it served as a good functional representation for most of the spectra analyzed, 
making it useful for comparison purposes. The reduced x 2 for the model-fits to 
the 4-channel data was, on average, higher than that for the 16-channel data; 
this is attributed to systematic errors. The uncertainties in the two parameters 
of interest, E ce nt and <7 g , were estimated jointly by asstuning a 68% confidence 
region in the shape of an ellipsoid, centered on the measured values. 5 

To address the issue of consistency of results using different datatypes we 
analyzed both 4-channel and 16-channel spectra integrated over the same time 
intervals for a small sample of bursts. In general, the best-fit parameters agreed 
within their uncertainties. The amplitude of the 16- channel data was usually 
slightly larger than that of the 4-channel data, which we attribute to slightly 
over- correcting the 16-channel data for dead-time effects. 

RESULTS 

Figure la shows a plot of o g vs £ cent for the set of 29 short bursts. Fig- 
ure lb shows the same information for the 55 short spikes within the 13 long- 
duration bursts. Though it is difficult to draw conclusions from such crowded 
plots, it does seem clear that there is a large concentration of short spikes within 
long-duration bursts at low values of E ce nt- The frequency distributions of E cent 
and <7g are given in Figures 2a and 2b, respectively. The distributions for the 
short- duration bursts are shown as solid lines, while the long-duration bursts 
are represented by dotted lines. Assuming normally distributed errors, the un- 
certainties for each bin are simply the square-root of the number of events in 
that bin. 
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Fig. 1 (a and b). la shows a plot of the best-fit model parameters, along with 
their uncertainties, for the short- duration bursts, lb shows the same information 
for the short spikes within the long-duration bursts. 


Figure 2a shows that the short spikes within the long-duration bursts and 
the short-duration bursts sample the same range of parameter space of E ceni . 
However, their distributions are clearly different; i.e., the short spikes in long- 
duration bursts tend to concentrate at low values of E cent while the values for 
the short bursts seem to be distributed uniformly over the entire range. This 
is in contrast to the result of comparing the distributions of <r» for each type of 
event, shown in Figure 2b. It indicates that the short bursts ao not sample the 
same range of available parameter space as the short spikes, and moreover, that 
the range of values of cr* for the short bursts is a subset of those sampled by the 
short spikes. Further, there is a hint of a bimodality in the distribution of <r g in 
Figure 2b for the short spikes, thus suggesting the interesting possibility of two 
classes of short spikes within long-duration bursts, based on spectral parameters 
alone. At this point, no work has been done to identify which long-duration 
bursts would makeup the two populations shown in Figure 2b, nor has any 



Fig. 2 (a and b). 2a and 2b show the frequency distributions of E mtx and <r g , 
respectively. 





attempt been made to correlate <r g with any other properties (peak intensity, 
duration, sky distribution, etc.) of the spikes or their “parent” bursts. 

It was found that the three freely varying parameters of the function used 
to model the photon spectrum are significantly correlated. Therefore it is not 
appropriate to independently use any one of these parameters to characterize a 
spectrum. To reduce the impact of this correlation, the energy at which the v F„ 
spectrum is a maximum, E m&X) was calculated analytically for each burst. This 
is done by multiplying Equation 1 by a factor of E 2 and setting the derivative 
of the result with respect to E equal to zero. Solving for E maz gives 

2,7 2 

logio(-^max) = log 10 (E l cen t) + : T ( 3 ) 

logio(e) 


•Em ax is then used to characterize each spectrum, with the hope that the effect 
of the correlation between the parameters has been minimized. 

The resulting distributions of E max for the two types of bursts are similar 
and are shown in Figure 3. The mean value of £ max for the short bursts is 
394 keV while for the short spikes within the long bursts it is 292 keV. This 
is consistent with the previously- known result that the short bursts tend to be 
harder than the long bursts. 2 The uncertainties for each bin in Figure 3 may be 
approximated as the square-root of the number of events in that bin, assuming 
that the errors in £ mix are distributed normally. The uncertainty in £ max was 
calculated in a standard fashion. 6 Assuming that the errors in E ceat and cr g axe 
small, the uncertainty in E ma , x was found using 
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where SE cent and 8cr* axe the uncertainties in E cent and <r g , respectively, and the 
partial derivatives ot E m&x axe calculated using Equation 3. The relative errors 
5-30% 


were usually in the D-< 3 U 7 o range 
The uncertainties in E max were of- 
ten at, or smaller than, the energy 
resolution of the data. Therefore we J2 
do not expect the uncertainties in S 
-fi'max to significantly change the re- u 
suits shown in Figure 3. 3 

With a few simple assump- 
tions, the results shown in Figure 3 
may provide some physical insight 
into the spatial distribution of burst 
sources. Let us assume that burst 
sources are at cosmological distances 
and that the mechanism producing 
the short-duration bursts is the same 
mechanism that produces the short 
spikes within the long-duration 
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bursts. By using the mean peak energy of the v F* spectrum for each distribution 
as a characteristic energy, we have calculated the cosmological redshift of the 
long-duration bursts zi relative to the redshift of the short-duration bursts zs . 
We find that under these assumptions zi, = 1.352s- However, a cosmological 
interpretation is not required to explain Figure 3; the fact that short bursts tend 
to be harder than long bursts, convolved with any tendency for hard-to-soft 
spectral evolution, would give the same results. 

A cautious approach must be taken in interpreting Figures 2 and 3 due 
to the low statistical significance of the data. The uncertainties in the best-fit 
parameters may not be normally distributed; thus, the errors on the frequency 
distributions should be calculated directly from the best-fit parameter uncer- 
tainties. It is clear that the uncertainties by themselves will not be able to wash 
out the strongly peaked distribution of E cent for the short spikes within long- 
duration bursts shown in Figure 2a. However, the bimodal distribution of er g 
must not be taken too seriously until the uncertainties axe more closely exam- 
ined. Also, it should be pointed out that Equation 4, used to calculate 6E ma , x , 
breaks down when the uncertainties 6E ceat or Scr s are large because second-order 
terms in the Taylor expansion have been neglected. 

FUTURE WORK 

There are a number of ways to improve the procedure used in this study. 
An attempt to choose a larger sample of bursts will be made, allowing more 
confidence in the interpretation of any interesting results. A function must be 
chosen for the photon spectrum model in which the freely varying parameters 
are less correlated. The limitation here is that the only BATSE data types 
available to study the short-bursts often consist of 4-channel energy spectra; 
thus the model must contain no more than 3 free parameters to ensure at least 
one degree-of-freedom in the fitting procedure. It may be possible to use a more 
complicated model with one parameter fixed to a value that is fairly represen- 
tative of most of the bursts being studied. Additionally, the assumption of a 
symmetrical ellipse for the shape of the 68% confidence region when determin- 
ing the uncertainties in the best-fit parameters may not be accurate. Software 
to examine these regions for each fit is in place and will be utilized in the fu- 
ture. Finally, the uncertainties in the best-fit parameters should be propagated 
through the frequency distributions in a more rigorous way. This will indicate 
how serious one can take features such as the bimodality in Figure 2b. 
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ABSTRACT 

One of the primary objectives of BATSE is to determine the direction to 
each observed gamma- ray burst, allowing the angular distribution of the burst 
sources to be studied. Most of the BATSE burst data types are not suitable 
for this purpose, since they are produced by combining rates from several detec- 
tors (details of the BATSE experiment and the various datatypes can be found 
elsewhere 1,2 ). After the loss of the CGRO tape recorders in March, 1992, the 
availability of data with rates from individual detectors has often been limited 
by telemetry gaps. Consequently, many burst directions are derived using the 
MAXBC data, which are almost always available, but have significant limita- 
tions. MAXBC data contain the maximum background-subracted count rate 
(counts per 64 ms) in the energy range of ~ 50-300 keV, averaged over one 
second, for each of the eight Large Area Detectors. It is important to evaluate 
the locations computed using MAXBC data because these locations will appear 
in the BATSE 2B Catalog when no other data are available. We report on the 
comparison of burst locations computed using MAXBC data with locations de- 
rived from the Third Interplanetary Network 3 , hereafter known as IPN 3 . We 
report on the comparison of MAXBC locations of hard solar flares with the 
known position of the Sun. We also report on the comparison of burst locations 
computed using MAXBC data with those locations given in the first BATSE 
catalog, for which the burst data are essentially complete. 4 

METHODOLOGY 

The 11 gamma- ray bursts used in these comparisons were selected because 
each has been well-localized by the IPN 3 . The 22 solar flares were chosen because 
their spectra are harder than an average flare, thus closer to a typical burst 
spectrum. Moreover, the flares extend the intensity range of our sample of well- 
localized events further down than possible with only the IPN 3 bursts. These 33 
events are the same events used in the analysis shown in Figure 2 of the BATSE 
IB Catalog. 

Three locations were calculated for each event. Details concerning the 
method of calculating locations, and the systematic and statistical errors in- 
volved, will not be discussed here but can be found elsewhere. 5 ’ 6 One of the 
three locations was determined using a power-law spectral index of 2.0. This 
value is used by the BATSE team when calculating locations of events for which 
there is no other suitable LAD data. Additionally, a value of 2.0 is used by 
the BATSE team when determining the preliminary location for a strong burst, 
with the intention of quickly alerting investigators observing at other wave- 
lengths. The second location computed for each event employs a spectral index 



( 1 ) 


obtained from a hardness ratio HR look-up table, with 

10000 xCh2 
HR ~ Ch2+ ChZ 

where Ck2 and ChS are the background-subtracted count rates in the LAD 
discriminator channel 2 (spanning ~ 50-100 keV) and discriminator channel 3 
(spanning ~ 100-300 keV) respectively. This method is used by the BATSE 
team when the preferred 16-channel CONTINUOUS data and 4-channel DIS- 
CLA data are not available due to telemetry gaps. The third location computed 
for each event utilizes the best-fit spectral index used in calculating the locations 
in the BATSE IB Catalog. 

The offsets (in degrees) of each location from the well-known locations 
were calculated. The offsets of each location from the published IB Catalog 
location (for which the burst data are essentially complete) are also calculated. 
The values used for the IB Catalog solar flare locations are the same values 
used to make Figure 2 of the IB Catalog. In addition to locations, the signal- 
to-noise ratio in the 2 brightest detectors was calculated for each event. If one 
trys to compare Figure 2 of the 2B Catalog with the results presented here, it 
must be noted that the value of the signal-to-noise ratio for each event will not 
necessarily be the same. In the IB Catalog we were able to use various data 
types, most of which provided better energy resolution than MAXBC, as well as 
more flexibility in choosing the burst time interval (recall that MAXBC are the 
background- subtracted rates only at the peak time of the burst in the 50-300 
keV energy range). 


RESULTS 

Table I. provides a brief summary of the various location comparisons. 
The two samples being computed are given in the first two columns. The third 
column specifies the method used to obtain the power-law sepctral indices. The 
fourth column gives the mean total (systematic + statistical) error (a ) , where 

( 2 ) 

1=1 

and <7j is calculated as the difference (in degrees) between the locations from 
each of the two samples, for the event. N is the total number of events 
in the sample. The fifth column informs the reader which figure contains the 
corresponding comparison data. 

Figure 1 shows the offset of the MAXBC locations from the well-localized 
locations plotted as a function of the signal-to-noise ratio in the two brightest 
detectors. The distinguishing characteristic for each of the three plots is the 
method used to obtain the power-law spectral indices. The dashed line in each 
figure is a representation of the statistical errors as calculated in the BATSE IB 
Catalog. 

Figure 2 shows the offset of the MAXBC locations from the locations 
given in the IB Catalog plotted as a function of the signal-to-noise ratio in the 
brightest detectors. Again, the distinguishing factor between these three plots 
is the method used to obtain the power-law sepctral indices. The dashed line 
again represents the BATSE IB Catalog statistical errors. 
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FUTURE WORK 


These results, combined with a more detailed analysis of the seperate 
systematic and statistical errors, will be included in the BATSE 2B Catalog. 


j Table I. Analysis of MAXBC Locations j 
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ABSTRACT 

A study of the continuum spectra of gamma-ray bursts (GRBs) and solar 
flaxes (SFs) is performed to investigate the photon energy flux of these events. 
BATSE 16 channel Large Area Detector (LAD) data are fit using a forward- 
folding model-fitting technique to produce photon spectra. The integral of the 
fit of each spectrum is used to compute the energy flux over consecutive intervals 
of the event, resulting in a time history of the power output (ergs /sec). These 
are compared to search for similarities among the emission mechanisms in GRBs 
and SFs. Comparison with the solar flare energy emission may yield information 
about the gamma-ray burst production site. 

DATA AND ANALYSIS 

The Burst and Transient Source Experiment (BATSE) consists of eight 
independent detector modules mounted on the corners of the Compton Gamma- 
Ray Observatory (CGRO), and is capable of nearly full-sky observations. This 
detector system, explained in detail elsewhere 1,2 , provides continuous monitoring 
of the x-ray and gamma-ray sky. Each of the detector modules consists of two 
Nal(T^) scintillation detectors: a Large Area Detector (LAD), optimized for 
temporal resolution, and a Spectroscopy Detector (SD), optimized for energy 
resolution. 

We use the BATSE Continuous (CONT) and Medium Energy Resolu- 
tion (MER) datatypes from the LADs. The CONT datatype, received in every 
telemetry packet, provides 2.048 second resolution in 16 energy channels span- 
ning the range ~20-1800 keV; the MER datatype, accumulated only in burst 
mode and beginning ~30 ms after a trigger occurs, consists of 16 energy chan- 
nels providing 16 ms resolution for ~32.8 seconds followed by 64 ms resolution 
for ~130 additional seconds. Thus MER data does not provide any pre-trigger 
information. Events which had significant emission before the trigger were, in 
general, analyzed with CONT data; events with fine time structure were ana- 
lyzed with MER data to enable resolution of shorter spikes within longer events. 

The data were deconvolved using a standard forward-folding model-fitting 
technique to produce photon spectra. Each event was divided into time intervals 
of equal duration and a photon spectrum was produced for each interval. The 
spectra were integrated over the fitted energy range (~20-1800 keV for GRBs 
and ~30-225 keV for SFs) to obtain the total photon and energy output during 
each interval, resulting in a time history of the photon and energy flux for that 
event. The flux histories of gamma-ray bursts and solar flares were compared to 
search for similarities among the two types of events. We find that the solar flare 
photon and energy flux histories are of almost identical shape, indicating little or 
no spectral evolution for these events. The gamma-ray bursts sometimes show 
significant differences between the photon flux and and the energy flux. This 
gives some information regarding the emission mechanism of the burst sources, 



suggesting that the burst emission environment is not the same as that of solar 
flares. 

The gamma-ray burst data were fit with a normal distribution in logio 
energy ( E) with variable (linear function of energy) width over the —20-1800 
keV range: 


dN 

dE 




'loglo(-E) ~ logio (Ecent) 
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( 1 ) 


where 


and 


logio (FWHM(E)) 
2.35482 


logio (FWHM(E)) — cr(Ecent) + Ac(£^cent) X (logio(E) — logio(Ecent)) • 


There are four varying parameters in this model: A denotes the amplitude, E ce n t 
the energy of the centroid, cr( JE7 ce nt) the full width at half maximum ( FWHM ) 
of the Gaussian, and Acr(J5J cent ) the change in the full width with energy. This 
function fits the wide range of observed gamma-ray burst spectra well, keeping 
in mind that we are not strictly concerned with the model that is used, but only 
that the function represents the data reasonably well. Therefore this functional 
form is adequate to ensure that the integral is descriptive of the actual energy 
output. The solar flare data were fit with a simple power law over the —30-225 
keV range: 



with Eo fixed at 100 keV. A is the amplitude of the function, and a is the 
power-law index (generally —3-6 for the solar data examined here). Limiting 
the energy range over which the solar flare data were fit enabled a straight power 
law to sufficiently describe the spectra of these events. The power law fits to 
the solar data sometimes exhibited large x 2 during the most intense portions of 
the events due to the excellent counting statistics at these times. 

RESULTS 

The analysis was performed for 15 solar flares and 11 gamma-ray bursts 
observed by BATSE. Although the sample sizes were limited, we attempted to 
select events which differed significantly in their time histories. Eight of the 15 
solar flares were visible above —100 keV, the remainder showing emission only 
below -100 keV. Seven of the 11 gamma- ray bursts showed significant emission 
above —300 keV. 

For emission produced by a single mechanism, the absence of extreme 
spectral evolution implies that the energy flux history closely mimics the pho- 
ton flux history. The solar flares examined exhibited this property, supporting 
observations made with the Gamma-Ray Spectrometer (GRS ) 3 on the Solar 
Maximum Mission (SMM), which showed that the hard x-ray (—10-140 keV) 
and 7 -ray (—10-100 MeV) fluxes peak within —one second of each other. It has 
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also been shown that the radio flux peaks near this time as well (see reference 4 
and references therein). Figure 1 shows two typical examples of the flares which 
were analyzed, demonstrating the excellent correlation between photon flux and 
energy flux. It should be noted that it is possible for the two varying parameters 
of the power law model used to fit the solar data to change in such a manner 
that the energy evolution follows the photon evolution. 

Several of the gamma-ray burst profiles appear to exhibit a more pro- 
nounced rise in the energy flux followed by a steady decay even though the 
photon flux may be increasing (see Figure 2aL The burst shown in Figure 2a 
reaches its maximum in energy at the start of the event, while the photon flux 
peaks at approximately 20 seconds later. The smooth gamma-ray burst shown 
in Figure 2b reveals particularly intriguing features in its energy output evolu- 
tion. The energy reaches maximum quickly, then remains relatively constant for 
several seconds, implying sustained energy injection or lack of energy dissipation 
in the burst environment. There is also a striking exponential decay in energy 
at the end of the event, although the count rate history also shows feature. 

CONCLUSIONS 

Although the photon and energy fluxes are strongly correlated, examina- 
tion of the energy output of solar flares and gamma-ray bursts suggests that the 
emission environments of the two types of events are, in general, not identical. 
Several of the gamma-ray bursts reveal intriguing differences between the pho- 
ton output and the energy output. Observations with the BATSE Spectroscopy 
Detectors in conjunction with those of the LADs support the hypothesis of two 
separate photon production mechanisms in solar flares: thermal (^10 keV) and 
thick-target (i£10 keV) Bremsstrahlung. The diversity in the temporal profiles 
of GRBs with energy suggests that these events are also the result of two or more 
photon production mechanisms, albeit mechanisms which appear to be distinct 
from those responsible for solar flares. Comparison of the solar and gamma-ray 
burst photon and energy flux histories support this suggestion. 
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ABSTRACT 

The continuum spectra of bursts are studied with particular emphasis placed 
on the search for spectral subcomponents. The spectra observed in individual 
bursts show significant evolution throughout their durations. This study addresses 
the degree to which hard and soft spectra are combined in bursts and the effect 
this has on the complexity of the observed spectra. In particular it is shown that 
the presence of broad cusps in spectra in the 40-100 keV range can be explained 
by the superposition of hard and soft spectra in the bursts. This interpretation is 
not meant to explain narrow spectral absorption lines in bursts. It is simply an 
alternate explanation for some of the broader features observed in the continuum 
spectra of GRBs. 

INTRODUCTION 

Gamma-ray bursts are known to exhibit a wide range on continuum spectra 1 . 
Spectral variation is not only seen between bursts but within bursts as well 2,3 . 
Bursts have been observed with narrow absorption and emission lines in them 
as well as with broader features of this type 1,2,4 ' 5,6 . This work shows that broad 
absorption features in the 40-100 keV range can be explained by the superposi- 
tion of hard and soft burst continuum spectra. These results are important for 
estimating the frequency with which absorption line features should be observed 
in gamma-ray burst continuum spectra. It also has implications for the temporal 
behavior of the energy release mechanisms operating in bursts. 

PROCEDURE AND DATA ANALYSIS 

The medium energy resolution (MER) data from the BATSE 7 instrument, 
consisting of 16 channel Large Area Detector (LAD) data in the 20-2000 keV 
energy range with up to 16 ms time resolution, were used for the bulk of this 
analysis. The LAD spectra were compared for consistency to the spectroscopy 
(SD) detector data and were found to be in agreement. 

The main spectral continuum model used here is a generalized form of the 
lognormal distribution. It is a gaussian of the log of the energy where the full 
width half maximum (FWHM) is a linear function of the log of the energy. This 



is a very flexible functional form that will fit spectra that do not have a change 
in the sign of their curvature in the interval being fit. Its mathematical form is 
presented in detail elsewhere 8 . 

The channel to energy conversion for the LAD 16 channel data has been 
optimized 9 using in-flight calibration data to remove systematic effects from the 
spectral analysis. 

Figure 1 shows two very different spectra observed in the same burst fit with 
the model described above. Figure 1A shows the continuum spectrum for the first 
1.2 s of BATSE trigger 907 (IB 911007) while figure IB shows a 1.2 s spectrum 
117.5 s into the burst. The superposition of these two spectra would produce a 
cusp or broad absorption feature in the resulting spectrum. 





Energy (keV) Energy (keV) 


Figure 1: Spectral fits for two separate intervals 
of trigger 907 (IB 911007) 


The obvious next step is to inspect the time profiles of the bursts to find 
intervals where the onset of hard Figure 1A type emission is preceded by softer 
emission and see whether significant cusps are visible in the spectra. Figure 2 
shows the time profiles of the counts spectra of trigger 1025 (IB 911109) in four 
broad energy channels. These are the LAD discriminator data. In this burst 
there is an episode of hard emission 2s after the burst trigger lasting 0.3 s, as is 
evident in the E > 300 keV time profile. The 20-50 keV time profile shows a fairly 
gradual rise in counts from the trigger time until the hard outburst. Furthermore 
no feature corresponding to the hard episode is visible in the 20-50 keV range. 

In order to test for the significance of two spectral component superposition 
in any of the spectra comprising this burst the MER data were summed into 
0.2s intervals over the duration of the event. These spectra were then fit with 
the generalized lognormal described above for each interval. Then the Sunyaev- 
Titarchuk Comptonization model was added to the generalized lognormal and 
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Figures 3 A and 3B: Spectral fits for two separate intervals 
of trigger 1025 (IB 911109) 


i mi 






the fitting was repeated. The F-test was then performed on each interval to 
see if the addition of the Sunyaev-Titarchuk Comptonization model produced a 
significant improvement in the fit. The improvement was significant at the 1% level 
(i.e. a 1% chance that the improvement was incidental) for the 2. 0-2. 2s interval 
and at the 5% level for the two intervals preceding it. All the other intervals 
showed no significant improvement even though many of these intervals had source 
count rates comparable to the intervals showing two component structure. The 
16 channel, 2.048 s resolution LAD data for each separate detector showed the 
significance of the added continuum component as well. The SD detector data 
were consistent with the LAD data although not as statistically significant. 

Figure 3 A shows a spectrum from the 1.2- 1.4 s interval of trigger 1025 fit with 
the generalized lognormal acceptable at the 10% level (i.e. the model will produce 
a data set like this at least 10% of the time). Addition of the Sunyaev-Titarchuk 
model yields no significant improvement. Figure 3B shows the 1. 9-2.1 s interval 
of trigger 1025 fit with the generalized lognormal and the Sunyaev-Titarchuk 
model acceptable at the 10% level. The generalized lognormal fit alone is only 
acceptable at the 5.0 x 10~ 5 % level. The F-test shows that the addition of the 
Sunyaev-Titarchuk model is a significant improvement for this interval at the 1% 
level. 

Another burst that shows some broad cusp-like structure in the first second 
after reaching trigger threshold is 1085 (IB 911118). Figure 4A shows the photon 
spectrum for the first 0.0-0.25 s interval of the burst fit to the generalized lognor- 
mal. Figure 4B shows the counts spectrum for the same fit to demonstrate the 
presence of the structure in counts space. Although the generalized lognormal fit 
in the 20 to 300 keV range cannot be ruled out since it is acceptable at the 8% 
level, the F-test shows that the addition of the Sunyaev-Titarchuk model is a sig- 
nificant improvement for this interval at considerably greater than 1% level. The 
pattern of evolution of this burst over the first second since trigger is one where 
a hard component grows in intensity until it buries the softer emission obscuring 
the cusp-like feature present in the first quarter second. 

DISCUSSION 

This is not the first time that multiple spectral components have been pro- 
posed for bursts 10 or the first time that the superposition of spectral components 
has been proposed as an explanation for some of the broad, absorption line-like 
features observed in bursts 11 . The BATSE LAD data provide us with an opportu- 
nity to study the evolution of spectral components and absorption lines with good 
statistics and see which explanation fits the observations most convincingly. It 
seems unlikely that the superposition of spectral components could explain narrow 
absorption features that have been observed with other instruments 5 , however it 
could explain some of the more complicated looking spectra found in the Konus 
Catalog 1 . Determining how many of the previous observations of lines are clearly 
narrow absorption lines vs. broad cusps is important for determining how likely it 



tinuum spectra. Studying the comparative evolution of hard and soft components 
in bursts may yield valuable insights into the physical mechanisms that produce 
the observed spectra. For instance the time profile of trigger 1025 clearly shows 
the onset of hard emission after the burst process has been initiated. This shows 
that the bursts do not evolve uniformly from hard to soft and indicates that hard 
emission process can trigger repeatedly during a burst’s output. Studying the 
development of the bursts’ spectra should help map out the characteristics of this 
as yet unexplained physical phenomenon. 




Figures 4A and 4B: Spectral fits for the same interval of 
trigger 1085 (IB 911118) in units of photons and counts. 


REFERENCES 

1. Mazets et al, 1981, Astophys. Space Sci. 80:3-143 

2. Mazets et al, 1982, Astophys. Space Sci. 82:261-282 

3. Mitrofanov et al, 1984, Sov. Astron. AJ 28:547-49 

4. Teegarden et al, 1980, Ap J. Lett. 236:L67-70 

5. Murakami et al., 1988, Nature 335:234-35 

6. Golenetskii et al 1986,Astrophys. Space Sci. 124:243-78 

7. Fishman et al, 1989 in proceedings of the GRO Science Workshop, 
ed. N. Johnson, pp. (2)39-2(50) 

8. Mallozzi et al, these proceedings. 

9. Pendleton et al, proceedings of the second compton symposium. Sept 20-22, 
1993, U. Maryland, in press. 

10. Mitrofanov et al, 1992 in Gamma-Ray Bursts: Observations, Analyses, and 
Theories, Ed. C. Ho, R. Epstein, and E. Fenimore, pp. 209-216 

11. Fenimore et al, 1982 in Gamma Ray Transients and Related Astrophysical 
Phenomena, Ed. R. Lingenfelter, H. Hudson, and D. Worrall 





0 ) ( 2 ) 


CN22 < 

TURNER J/PUBLICATION 
MARSHALL SPACE FLIGHT CENTER 
HUNTSVILLE AL « 


DELETIONS OR CHANGES 544-4494 
RETURN ADDRESS CN22D 000002444 



